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Abstract 

A high performance distihation system to remove krypton from xenon was con- 
structed, and a purity level of Kr/Xe = ~ 3 x 10~^^ was achieved. This develop- 
ment is crucial in facilitating high sensitivity low background experiments such as 
the search for dark matter in the universe. 
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1 Introduction 

Liquid xenon is one of the most attractive materials for use in detectors for 
astroparticlc and particle physics [1,2,3,4,5]. As a scintillator it has a large 
light yield, comparable to that of Nal(Tl), and can be used for detectors 
with low energy thresholds and high energy resolution. Because of the high 
atomic number of xenon {Z = 54) and its high density in liquid form (~ 3 
g/cm^), it contributes to the reduction of environmental backgrounds, such as 
7-rays and /3-rays from uranium and thorium contamination, by self-shielding. 
Another big advantage of liquid xenon is that xenon does not have long- 
lived radioactive isotopes, and thus experiments on rare phenomena (such as 
dark matter searches and double beta decay experiments) may be carried out 
shortly after moving the xenon underground from the surface. 

One drawback of liquid xenon is the possibility of contamination with radioac- 
tive ^^Kr. Xenon is usually produced from air. The concentration of xenon in 
air is ~ 10~^ mol/mol, while the concentration of krypton is ~ 10~^ mol/mol. 
In the commercial production of xenon, krypton is removed by distillation 
or adsorption. But, the commercially available xenon contains 10~^ ~ 10~^ 
mol/mol of krypton, because the removal starts from the Kr/Xe ratio of about 
10 and such purity is enough for most of the applications of xenon. ^^Kr is a ra- 
dioactive nucleus which decays into rubidium-85 with a half-life of 10.76 years, 
and emits /9-rays with a maximum energy of 687 keV and a 99.57% branching 
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Fig. 1. Background event rate of ^^Kr for a Kr/Xe ratio of 10~^ mol/mol, compared 

with dark matter signal rate, assuming a cross section of 10~'^ pb and a WIMP mass 
of 100 GeV/c^. A quenching factor of 0.2 is assumed for the dark matter signaL 



ratio. (The remaining 0.43% represent a /9-emission with a maximum energy 
of 173 keV followed by a 514 keV 7-ray emission.) Large quantities of ^^Kr 
are produced artificially in nuclear fission. This constitutes the main source of 
^^Kr in air. The concentration of ^^Kr in air is measured to be ~ 1 Bq/m^[6,7], 
which corresponds to *^Kr/Kr = ~ 10^^^. Assuming a Kr/Xe ratio of 10"''' 
mol/mol, the background event rate of ^^Kr is as shown in Fig. 1. (^^Kr/Kr of 
1.2 X 10~^^ is used here.) The expected dark matter event rate is also plotted, 
assuming a cross section of 10"^ pb and a WIMP mass of 100 GeV/c^. A 
quenching factor of 0.2 is assumed for the dark matter signal. 

Because next-generation dark matter detectors aim at a sensitivity of 10^^ pb 
for the dark matter-proton cross section, a background rate of 10^"^ events 
/keV/day/kg is required. This corresponds to a Kr/Xe ratio of less than ~ 
10-^2 mol/mol. 

Possible methods to remove krypton from xenon are distillation and adsorp- 
tion. They are commonly used industrial processes but systems which re- 
duce Kr enough to meet our requirements did not exist. A development of 
adsorption-based chromatography to achieve the required Kr concentration is 
reported in ref.[8]. In this paper, we describe the development of a distillation 
system to reduce krypton to a level of 10^^^. In section 2, the design princi- 
ple is discussed, while in section 3, the setup and operation of the distillation 
system are described. In section 4, we describe a technique for measuring low 
levels of krypton in xenon gas, and discuss the results. 
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Fig. 2. (a) Illustration of the McCabe-Thiele (M-T) method. The various elements 
of the figure are explained in the text, (b) Calculation of the number of theoretical 
cells using the M-T method. The horizontal and vertical axes represent the krypton 
concentration in the liquid and gas phases, respectively. The thick solid curve is 
the equilibrium curve, and the thick solid lines are the condensation line and the 
collection line (details explained in the text). 

2 Design principle 

The boiling point of liquid krypton is 120 K at 1 atmosphere, while that 
of xenon is 165 K. This means that, in principle, separation of krypton and 
xenon can be performed by a distillation method. For the basic design of the 
distillation system, we followed the McCabe-Thiele (M-T) method [9]. It is the 
standard method of designing a distillation system. At first, the concept of 
the M-T method is reviewed, and then the application of the method for our 
purpose is described. 

The principle of the M-T method is shown in Fig. 2. The main element in the 
distillation system is a tower in which gas-liquid equilibrium is maintained. A 
vessel called a "reboiler", at the bottom of the tower, boils the liquid using 
a heater. In order to maintain a constant temperature profile in the tower, 
a condenser is placed at its top. The supplied xenon gas is cooled down to 
near boiling point, and then supplied to the feed point in the tower at a flow 
rate F, as shown in the figure. The processed xenon, which contains a lower 
concentration of krypton than the original xenon, is obtained from the reboiler 
(flow rate W), and xenon with a higher krypton concentration is obtained at 
the top of the tower (flow rate D). In the following, the Kr concentration of 
xenon at these points is expressed and xd, respectively. The heating 

power of the reboiler and the cooling power of the condenser at the top of the 
tower control the flow of xenon in the tower (flow rate L). The reflux ratio. 
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represented hj R = L/D. indicates the amount of xenon returned to the tower 
compared with the amount extracted with a high concentration of krypton. 

In the M-T method, the distillation tower is assumed to be a connected series 
of theoretical cells, with the gas-liquid equilibrium changing by one step in 
each cell. The number of theoretical cells is calculated for given boundary 
conditions of F, D, W, xp, xw, xd and R. 

For an equilibrium of liquid and gas phases in a mixture of elements (in this 
case, Kr and Xe), the partial pressures of the gas-phase elements (expressed 
as pKr and pxe, respectively) can be related to the fraction of each element in 
the liquid phase (expressed as xkv and xxe) by Raoult's law: 



PKr = PKr ■ XKr (1) 
PXe = Pxe ■ XXe (2) 

where Pkt and Pxe are the vapor pressures of the corresponding elements 
(2090 kPa and 201.4 kPa, respectively, at 178 K). Using Raoult's law, the 
concentrations of Kr in the gas phase {yKr) and the liquid phase {xkt) can be 
related by the equation 



VKt = —7 7T 3 

1 + {a - l)XKr 

where a is the ratio of vapor pressures of Kr and Xe, Pxr/Pxe (10.4 at 178 
K). Eq. 3 is called the "equilibrium curve", and it is shown by the thick 
solid curve in Fig. 2(b). For each theoretical cell, the gas- and liquid-phase 
Kr concentrations in the neighboring cells are related by conservation of mass 
flow, and may be expressed as: 



yn+l^^X. + ^X, (4) 

R' 1 ... 

where eq. 4 (5) is for the cells above (below) the feeding point, and n (m) is 
the number of theoretical cells. R' is (L -|- qF)/W, where q is the fraction of 
hquid with respect to the total feed material. Equations 4 and 5 represent the 
"condensation line" and the "collection hue" , respectively; these are shown as 
solid lines in Fig. 2(b). Figure 2(b) illustrates the method of estimating the 
number of required theoretical cells. 

We designed the distillation tower to fulfill the following requirements: 
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• The Kr concentration of the processed xenon should be three orders of 
magnitude smaller than that of the original xenon sample, i.e. xw = x 

Xp. 

• The collection efficiency of xenon should be 99%, i.e. W/F = 0.99 and D/F 
= 0.01. 

• The system should have a process speed of 0.6 kg Xe per hour, which allows 
100 kg xenon to be purified within one week. 

• The system should have a refiux ratio of R = 191, which means the heating 
power required for the reboiler is 14 W. 

• Xenon should be fed into the system in the gas phase, i.e. g = 0. 

A M-T diagram based on these requirements is shown in Fig. 3, in which the 
concentration of Kr in original X.e{xF) was assumed to be 3 x 10~^, the actual 
value of the processed xenon described in section 4. The M-T diagram shows 
that we need 6 stages of theoretical cells. 




Fig. 3. M-T diagram with system requirements. 



3 Setup and operation of the distillation system 

A key element of the distillation system is the packing column in the distilla- 
tion tower. We used "DX laboratory packing" (Sulzer Co.) with a diameter of 2 
cm. The catalog specifications of the packing included a liquid load (defined as 
the fiow per unit area per unit time) between 0.1 and 5 m^/rn^h. In our case, 
the liquid load is 1.2 /in^h for a refiux ratio R of 191. Gas loads, defined as 
Vyfp (where v and p are gas velocity and density, respectively), are 0.11 Paa 
and 0.23 Paa for below and above the feed point, respectively. The catalog of 
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the packing gives the height equivalent to a theoretical plate (HETP) for the 
distillation of a ethylbenzene/chlorobenzene mixture. The given HETP value 
is 2~8 cm for the gas load of 0.2-1.2 Pa2. However, since the HETP value 
might strongly depend on the liquid/gas loads and the type of liquid used, we 
conservatively enlarged the HETP value to 35 cm for the Kr/Xe distillation in 
the design, and the total length of the tower is set to 208 cm. The feed enters 
the tower 10.8 cm from the top. 

The full setup of the distillation system is shown in Fig. 4. Before entering the 

condenser 




Fig. 4. Flow chart of distillation system. 

distillation apparatus, the xenon flows through a getter SPRG-100H-00030X 
(Taiyo Toyo Sanso Co.), which is able to decrease concentrations of N2, CH4, 
O2, CO2, and CO to below 1 x 10"^° mol/mol, and H2 and H2O below 5 x 10"^° 
mol / mol. The xenon then flows through a heat exchanger which pre-cools the 
xenon. For this purpose, the purified xenon gas extracted from the distillation 
apparatus is used as a cooling medium. After the pre-cooling, the xenon is 
further cooled down to 192 K by a GM-type refrigerator (TZ-65; Taiyo Nippon 
Sanso Co.) which has a cooling power of 100 W at around 180 K, and the cooled 
xenon is fed into the distillation tower. The temperature of the condenser at 
the top of the distillation tower is kept at 178 K by an electrical heater. 
The reboiler is a cylindrical copper vessel with a 200 mm diameter. Three 
liquid level sensors are placed in the reboiler at 5 mm {LL level), 20 mm (L 
level) and 65 mm {H level), corresponding to 0.47 kg, 1.88 kg and 6.1 kg of 
xenon, respectively. The heater attached to the reboiler has a power of 22 W, 
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and runs with a 65 % duty cycle during operation (this corresponds to 14 W 
effective power) . The purified xenon is collected from the reboiler and goes to a 
collection bottle via the heat exchanger. Xenon with a high Kr concentration, 
which is called the "off-gas" , is collected from the top of the distillation tower. 
The purified and off-gas xenon were collected in stainless steel bottles which 
were cooled by liquid nitrogen. 

The apparatus was used from March 9 to 17, 2004, to process 100 kg of xenon. 
First, about 2.0 kg of xenon were liquefied and supplied to the distillation 
tower up to liquid level L in the reboiler. The heater in the reboiler and the 
condenser at the top of the distillation tower were switched on, and the system 
was maintained for 11 hours without collecting purified and off-gas xenon in 
order to establish the initial gas-liquid equilibrium in the tower. After this 
time, the fiow rate of the purified xenon was set to 0.6 kg/h (1.7 L/min gas 
fiow), and that of the off-gas xenon was set to 0.006 kg/h (17 mL/min). The 
flow rate of the xenon feed was set at around 0.6 kg/h (1.7 L/min) in the 
following manner: after the liquid level in the reboiler had reached L, the rate 
was set to 0.72 kg/h (2.0 L/min) for about four hours, and then to 0.48 kg/h 
(1.4 L/min) until the xenon in the reboiler returned to level L. When we 
found that the liquid level around L could be measured more precisely by the 
monitored pressure of the purified xenon, the fiow rate of the xenon feed was 
fine-tuned to balance the output fiow. Figure 5 shows the time variation for 
off-gas, input and purified xenon. Figure 6 shows the temperature variation as 
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Fig. 5. Time variation of flow rate of off-gas xenon (upper figure), and input (middle) 
and purified xenon (lower) . 

a function of time for the condenser in the distillation tower, the reboiler and 
the center position in the tower. The temperature was stable to within ±0.5 
K during operation. The pressure in the distillation tower was approximately 
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03/09 03/10 03/11 03/12 03/13 03/14 03/15 03/16 03/17 

date 

Fig. 6. Temperature variation of condenser, (TI22), reboiler (TIC24) and center 
position in the tower (TIC23). 

0.105 MPa (gauge) throughout the operation. 



4 Measurement of the krypton concentration 

In this section we describe the different methods used to measure the krypton 
concentrations in the input xenon, off-gas, and purified xenon. 

4-1 Krypton concentration in the xenon feed 

The Kr concentration in the original xenon sample was evaluated by a scintilla- 
tion measurement in a prototype detector for the XMASS experiment [1]. The 
original xenon sample (100 kg) was placed in a 30 L chamber, and scintillation 
light was detected using 54 low-background photomultiphers[10]. The back- 
ground spectra were measured in December 2003 before the distillation[ll] and 
in August 2004 after the distillation [12] as shown in Fig. 7. The background 
rate around 100-300 keV before the distillation is higher than that after the 
distillation by ~ 5 x 10"^ /kg/keV/day. This background level corresponds to 
a Kr concentration of (2 ~ 3) x 10~^ mol/mol, assuming a ^^Kr/Kr ratio of 
1.2x10-^^ 



4-2 Krypton concentration in the off-gas 

The Kr concentration of the off-gas xenon was measured on March 23, 2004. 
Xe and Kr were separated by gas chromatography (GC), and the ion count 
was measured by a photoionization detector (PID). The GC system for this 
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Fig. 7. The background spectra measured in December 2003 before the distilla- 
tion(sohd histogram) and in August 2004 after the distillation (dotted histogram). 
The dashed curve shows expected ^^Kr beta spectrum assuming a Kr concentration 
of 3xl0~^ mol/mol. 

experiment was a GC-263-50 (Hitachi Co.) with a separation column consisting 
of a SUS tube of 3 mm inner diameter and 2 m length, which was filled with 
MS-13X molecular sieves (30/60 mesh). In order to calibrate the system, a 
sample gas consisting of 1 x 10"^ Kr, 1% Xe and 99% helium (1 ppm Kr 
gas) was measured using this setup, and the count of Kr at the PID is shown 
in Fig. 8(a) for a gas sample of 5 mL. The integrated count rate for the 
Kr peak was 20,426. The off-gas measurement is shown in Fig. 8(b) for a 
sample of 5 mL. The large peak in front of the Kr peak is due to the N2 
component in the off-gas. As described in the previous section, the off-gas 
was collected by a cold trap method using a collection bottle immersed in 
liquid nitrogen. Unfortunately, there was a small leak at the endcap of the 
off-gas collection bottle. The Kr concentration in the off-gas was evaluated by 
subtracting the the tail of the N2 contribution as shown in the figure, and the 
resulting integral count rate of Kr in the off-gas was 5858. Normalization of 
this value based on the calibration gas gave a Kr concentration in the off-gas of 
(3. 3 ±1.0) X 10~^ mol/mol. The error estimate includes uncertainty introduced 
by the N2 subtraction. 



4-3 Krypton concentration in the purified xenon 



A highly sensitive counting method was needed to measure the Kr concentra- 
tion in the purified Xe, because the concentration was expected to be as low 
as 1 X 10~^^ mol/mol. The measurement described below was carried out in 
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(a) 1 ppm Kr 



(b) off-gas sample 




Retension time 



Retension time 



Fig. 8. Kr count by PID detector. The horizontal axis represents travel time through 
the GC column, (a) Count rate for calibration gas consisting of 1 x 10^^ mol/mol 
Kr, 1% Xe and 99% helium, (b) Count rate for the off-gas sample. 

September 2004. We first describe the principle of the counting method. 

The most sensitive method of gas component measurement which can be ap- 
plied to the measurement of Kr concentration is API-MS (atmospheric pres- 
sure ionization mass spectroscopy). In this method, a sample gas can be in- 
troduced at near atmospheric pressure, which enables large amounts of gas 
to be loaded and consequently allows highly sensitive measurements. The gas 
that carries the sample into the API-MS system (the carrier gas) must have a 
greater ionization energy than the target elements. The carrier gas mixed with 
the sample gas is ionized at the first stage (at this time, most of the ionized gas 
is carrier gas), and a charge exchange interaction in the next stage transfers 
charge from the carrier gas to the target elements. The ionized elements are 
then extracted and measured using a mass spectrometer. We used helium as a 
carrier gas, because the ionization energies of He and Kr are 24.6 eV and 14.0 
eV, respectively. Since the ionization energy of Xe (12.3 eV) is lower than that 
of Kr, Xe must be removed before introducing the sample into the API-MS. 
We used a GC technique for this purpose. In order to load a large amount of 
the sample gas, we used a concentration method, as detailed later. 

Figure 9 shows the setup for high-sensitivity Kr measurement. First, the sam- 
ple Xe gas was transferred to a 5 mL sample volume. The 5 mL sample was 
then transferred to the GC system using helium gas. The GC system used for 
this was GC-8APT (Shimadzu Co.) and the separation column consisted of 
a SUS tube of 3 mm inner diameter and 2 m length, which was filled with 



11 




API-MS 



Fig. 9. Setup of high-sensitivity measurement of Kr in xenon gas. 

MS-5A molecular sieves (30/60 mesh). When the Kr arrived at the output line 
of the GC, the line was connected to a metal filter (Nihon Seisen Co., pore 
size 0.003 fim) which was cooled by liquid nitrogen, and the Kr was trapped 
by the metal filter. The trapping efficiency was evaluated based on the results 
obtained for the cahbration gas, and was found to be as high as 80%. The 
processes described so far were repeated 100 times (concentration process), 
which enabled loading 500 mL (5 mL x 100 times) of the sample gas. After 
concentration, the metal filter was warmed to room temperature by removal 
of the liquid nitrogen vessel. The helium carrier gas was then used to transfer 
the Kr trapped in the metal filter to the API-MS. Another GC system was 
placed just in front of the API-MS to separate Kr to ensure its counting. The 
GC system used at this point had a separation column consisting of a SUS 
tube of 3 mm inner diameter and 2 m length which was filled with MS-13X 
molecular sieves (30/60 mesh). The API-MS apparatus was an API-200 (VG 
Gas Analysis Systems Ltd). The count rate of the processed and concentrated 
Kr sample is shown in Fig. 10(a). The M/Z value (the mass number of the 
element divided by the charge) was set to 84 in this measurement in order 
to count ^^Kr"*" ions. In order to check the stability of the API-MS system, a 
standard mixed gas composed of 1 x 10~^ Kr, 1% Xe and 99% He ("1 ppm 
Kr") was injected before the introduction of the sample gas. The six peaks 
before cycle number 1300 in Fig. 10(a) are due to the injection of the "1 ppm 
Kr" . The real Kr signal of the purified Xe was counted at around cycle number 
1470. Figure 11 shows a closer view of the region around the signal. The signal 
count was defined as the integral counts within I20 cycles of the peak, and was 
calculated to be 219 ± 14.8. The background count of the API-MS was eval- 
uated using the counts 50 cycles before and 50 cycles after the signal region; 
the value obtained was 0.70 ± 0.08 counts/cycle. After subtracting the back- 
ground contribution, the true Kr count of the purified Xe was calculated to be 
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Fig. 10. Kr count rate measured by API-MS for (a) purified Xe gas, (b) He carrier 
gas, and (c) 50 ppt Kr calibration gas. Kr counting for each gas was conducted 
at the cycle marked by the arrow. The high peaks appearing periodically before 
the real measurements represent injections of 1 ppm Kr calibration gas, which is 
composed of 1 x 10~^ mol/mol Kr, 1% Xe and 99% He, and are used as a reference 
for overall normalization. 
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Fig. 11. Kr count rate in purified Xe. Expanded view of Fig. 10 in the signal region. 

184.0 ± 15.8. In order to estimate the Kr contamination in the carrier gas, the 
concentration process was repeated only using the He carrier gas (i.e. without 
introducing the 5 mL xenon sample during each cycle). The API-MS count for 
the He carrier gas is shown in Fig. 10(b). Similarly, in order to evaluate the 
overall counting efficiency, which incorporates the collection efficiency of Kr 
and the counting efficiency of the API-MS system, a calibration gas contain- 
ing 5 X 10~^^ mol/mol Kr (50 ppt Kr) was processed using the same method; 
its count rate is shown in Fig. 10(c). The initial signal counts, background 
count rates and true signal counts (with background rate subtracted) for the 
purified Xe sample, the He carrier gas and the 50 ppt Kr gas are summarized 
in Table 1. The true signal is corrected for the variation of their respective "1 
ppm Kr" signals. Arbitrarily setting this normalization to 1.0 for the purified 
Xe, the corrections for the He carrier gas and the 50 ppt Kr become 0.92 and 
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0.88 respectively. 





Purified Xe 


He carrier gas 


50 ppt Kr 


Signal 


219±14.8 


409±22.2 


1235±35.1 


BG count/cycle 


0.70±0.08 


6.41±0.25 


1.16±0.12 


True signal(S) 


184.0±15.8 


95.6±25.8 


1341.7±40.7 



Table 1 

Signal counts, background count rates and true signal counts for purified Xe, He car- 
rier gas, and 50 ppt Kr gas. The true signal is {signal— BG/cyclex50cycles) / factor, 
where the factor is the relative counting rate of the "Ippm Kr" . 

The Kr concentration of the purified Xe was calculated using the following 
formula: 



5 X 10~^^ 

Kr concentration = {(Purified Xe) — (He carrier gas)} x — — 

'■^ ^ ^ ^ (50 ppt Kr) 

where the values in round parentheses are the true signal counts given in Table 
1. The value obtained for Kr concentration was (3.3±1.1) xlO~^^ Kr/Xe [mol/mol]. 



5 Conclusion 



A distillation system for removing Kr from Xe down to a concentration of 
~ 10~^^ Kr/Xe [mol/mol] was developed. The system was designed to reduce 
Kr concentration by three orders of magnitude with 99% Xe collection effi- 
ciency (i.e., the amount of Xe rejected is only 1%) and with a process speed 
of 0.6 kg Xe/h. This distillation system was used to purify 100 kg xenon gas 
containing 3 x 10~® Kr mol/mol. The off-gas was found to have a Kr concentra- 
tion of (3.3 ± 1.0) X 10"'' mol/mol, which is consistent with the transfer of the 
majority of Kr in the original Xe sample to the off-gas. The Kr concentration 
of the purified Xe was measured as (3.3± 1.1) x lO"-*^^ Kr/Xe [mol/mol] using a 
gas chromatography apparatus and a API mass spectrometer. These measure- 
ments confirmed that the distillation system can reduce the amount of Kr in 
xenon gas by up to three orders of magnitude. The achieved Kr concentration 
satisfy our requirements for the next generation of dark matter experiments. 
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